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Abstract. Over 300 extrasolar planets have been found since 1992, 
. showing that planetary systems are common and exhibit an outstand- 

(-H ' ing variety of characteristics. As the number of detections grows and 

as models of planet formation progress to account for the existence of 
I ■ these new worlds, statistical studies and confrontations of observation 

O ' with theory allow to progressively unravel the key processes underlying 

planet formation. In this chapter we review the dominant contribution 
^ ' of Doppler spectroscopy to the present discoveries and to our general 

understanding of planetary systems. We also emphasize the synergy 
of Doppler spectroscopy and transit photometry in characterizing the 
physical properties of transiting extrasolar planets. As we will see, 
^ , Doppler spectroscopy has not reached its limits yet and it will un- 

doubtly play a leading role in the detection and characterization of the 
first Earth-mass planets. 



1 Introduction 



o 

^' 

o 

, The question of the existence of other worlds has been present in human his- 
IL^ tory for millennia but it is only recently that scientific evidence has confirmed 
what many had anticipated: planets do exist and are common outside the So- 
lar System. The first robust detection of another planetary system came in 1992 
^ \ with the discovery of two terrestrial-mass planets orbiting the pulsar PSR 1257-1-12 
(IWolszczan fc Frail Il992h . Interestingly, this discovery did not receive all the at- 



X 



tention that could have been expected, probably because these two planets orbit 
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a "dead star" very different from the Sun and much less likely to host life in its 
vicinity. From this anthropocentric perspective, the major milestone in the search 
for extrasolar planets was the disco very in 1995 of 51 Peg b , the first extrasolar 
planet found to orbit a Sun- like star (jMavor fc Queloall995l ). Although of Jovian 
nature, 51 Peg b orbits at 0.052 AU from its parent star, a striking characteristic 
when compared to the giant planets in the Solar System, which all orbit beyond 
5 AU. This proximity has been a major surprise and a serious challenge to planet 
formation theories. 

Thirteen years after the discovery of 51 Pegb, over 300 extrasolar planets have 
been detected, including many fascinating system^. Five different techi iiques have 



contr ibuted to these discoveries: pulsar timin g (4 planets detec ted; e.g. IWolszczan 



1997 ) , Doppler spec troscopy (292 planets; e.g. Udrv et al]|2007 ) , photometri c tran 



sits ( 52 planets; e.g. Charbonneau et al]l2007 ). microlensing (8 planets; e.g. lGaudi 
20071). and direct imaging (4 objects with a mass possibly below 20 Mjup; e.g. 
Beuzit et al.l 120071 ). These observational techniques have considerably improved 



in recent years and Doppler spectroscopy, which has contributed the bulk of the 
discoveries so far, now allows the detection of planets with half the mass of Uranus. 

Likewise, since the discovery of 51 Pegb, planet formation theory has made 
rapid progress. To overcome the challenge posed by the unexpected properties of 
the newly found planets, two different formation models have been proposed: core 
accretion and disk instability. According to the core accretion model, dust grains 
coagulate to form planetesimals, which then accumulate to build up planetary 
cores. Planetary cores reaching the critical mass of 5-15 before the dissipa- 
tion of the gaseous protoplanetary disk subseq uently accrete significant am ounts 
of nebular gas and become giant planets (e.g. Lissauer fc StevensonI 2007h . The 
remaining solid cores merge through giant impacts to form terrestrial planets (e.g. 



Nagasawa et al.l 120071 ). In the alternative disk instability model, giant planet s 



form by direct fragmentation of the protoplanetary disk (e.g. lDurisen et aTl2007l ). 
Quantitative predictions based on the disk instability scenario are still sparse be- 
cause the simulations are computationally challenging and involve complex physics. 
In contrast, core accretion is now mature enough to allow for detailed calculations, 
and explicit comparisons with the observed population of extrasolar planets are 
possible. 

In this chapter we review the dominant contribution of Doppler spectroscopy 
to planet discoveries and to our general understanding of planetary systems. In 
Sect. [2] we describe the Doppler technique itself. Then, we present the observa- 
tional results and we discuss their interpretation within the current theoretical 
framework. In Sect. [3] we consider the results on giant planets for which an ex- 
tended statistics is available. In Sect. |4] we present recent results on low-mass 
planets of Neptune- and Earth-type, and we discuss the present limitations on 
Doppler precision. We end this review by describing in Sect. [5] the role played 
by Doppler spectroscopy in the characterization of transiting planets. All these 



^ See the Extrasolar Planets Encyclopedia, |http:/ /exoplanet.eu/ ' for an up-to-date census of 
the discoveries. 
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achievements and results are summarized in Sect. |6l where we also outline future 
perspectives. 

2 Doppler spectroscopy 

2.1 Principle 

Doppler spectroscopy is an indirect detection method which uses the starlight 
to measure the gravitational influence of a planet on its host star. Specifically, 
Doppler spectroscopy is based on the following key observations: 

f . In a planetary system, the star and the planet orbit their common barycenter 
according to Newton's law of gravitation and to the laws of motion. The two 
barycentric and the relative orbits have the same periods and eccentricities, 
but semimajor axes in the proportions a^, : ap : a = nip : nii, : (nii, -\- nip) , where 
is the mass of the parent star and mp the mass of the planet. The three 
orbits are coplanar and the orientations of the two barycentric orbits differ 
by 180° within that plane. 

2. According to the Doppler-Fizeau effect (hereafter simply the Doppler effect), 
the light emitted by a source approaching (receeding from) the observer is 
shifted towards shorter (longer) wavelengths. In its simplest relativist form, 
the Doppler formula writes 

. = ^ = -l±H^-l (2.1) 

where z is the so-called redshift, A and Aq are the observed and rest wave- 
lengths, respectively, K is the radial velocity, and V is the total velocity 
relative to the observer. In the In A space, the wavelength shift is indepen- 
dent of the rest wavelength and provides a direct measurement of the relative 
velocity between the source and the observer. 

3. The visible portion of the spectra of F, G, K and M dwarfs contains a mul- 
titude of metal absorption lines (Fig.[T|). These lines constitute a convenient 
benchmark to measure wavelength shifts through the Doppler effect. 

Putting these three observations together, Doppler spectroscopy consists in 
monitoring potential Doppler shifts in the spectra of stars with numerous absorp- 
tion lines. In the absence of other phenomena susceptible of shifting the stellar 
lines or of modifying their profile (Sect. I2.4.2p . measured Doppler shifts are con- 
verted into radial (line-of-sight) velocity variations and interpreted as the reflex 
motion of the star due to an orbiting companion (planet, brown dwarf, or star). 
Since the two barycentric and the relative orbits are closely related, measuring the 
reflex motion of the parent star along the line of sight gives access to some of the 
orbital parameters and characteristics of the companion. 
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Fig. 1. Blue part of an ELODIE spectrum of the planet-host star 51 Peg illustrating the 
multitude of absorption lines present in the visible spectra of solar-type stars. 



2.2 Planetary orbits and planet characteristics from radial velocities 

A planet in a Keplerian orbit induces on its parent star a perturbation of the form 

Vrit) = K[cos{iy{t)+uj) + ecos{uj)]+j (2.2) 
where K is the velocity semiamplitude 

27ra^ sin i 

K= — 2.3 

P(l-e2)i/2 

oj is the longitude of periastron, and 7 is the systemic velocity (velocity of the 
barycenter). Since the true anomaly, ^{t), depends on the orbital period (P), 
eccentricity (e) and time of passage at periastron (Tq), fitting a radial- velocity 
time series with the Keplerian model described above yields six parameters: K, e, 
w, To, P, and 7 (Fig.©. 

For a planetary system, the velocity semiamplitude is related to the masses of 
the two components through the so-called mass function 

P^ = ^K^l-e^r (2.4) 
(m-Vf -|- rrip)^ 27rG 
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Fig. 2. Illustration of the influence of each of the six parameters determining the Ke- 
plerian model of Eq. 12.21 The baseline model has P — W days, e = 0, if = 50 ms~^, 
w = 0°, To = 2.5 days, and j — —30 ms~^. These orbital elements are varied in turn, 
one on each panel. 
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which reduces to the expression of the planet minimum mass 

Kml'\l~e^f^ (2.5) 

under the assumption that nip <C (G is the universal gravitational constant). 
Introducing the same approximation into Kepler's third law yields an expression 
for the semimajor axis of the relative orbit 

a^apc^(^^y ml"p^/^ (2.6) 

Therefore, fitting radial-velocity data with a Keplerian model to account for the 
presence of a planetary companion gives four of the six orbital elements of the 
relative orbit (the longitude of the ascending node, fi, and the orbital inclination, 
i, remain unknown). When the mass of the central star can be estimated by other 
means, the Keplerian fit additionally yields a lower limit on the planetary mass 
and the semimajor axis of the relative orbit. Although the true mass of the planet 
can be significantly different from the minimum mass, the two values agree within 
a factor of 2 (rrip < 2 mp sin i) in 87% of cases. 

For a planet in a circular orbit around a solar-mass star, Eg. 12.51 simplifies to 

X[ms-i] ~ 28.57 mp sin i[Mjup] [yr] (2.7) 

Since the probability of detecting a planetary signal depends essentially on the 
value of the velocity semiamplitude, this expression indicates that Doppler mea- 
surements favor the detection of planetary systems with massive and short-period 
planets. Applied to the Solar System, Eq. 12.71 shows that Jupiter induces on the 
Sun a radial-velocity perturbation with a semiamplitude of 12.5 ms^^, while the 
Earth induces a perturbation with a semiamplitude of 9 cms~^. As we will discuss 
now, obtaining radial velocities at this level of precision is nontrivial. 



P 
2^ 



2.3 Instrumentation and techniques for high-precision Doppler measurements 



A measurement error of 1 ms~^ corresponds to AA/Aq = 3 x 10~^. On the detec- 
tor of a typical spectrograph (i?~ 100,000), this wavelength shift translates into 
a linear displacement of about 3x 10~^ resolution element (~10~^ line width). 
Although attaining this precision is feasible by averaging over many spectral lines, 
experience shows that unless one takes special measures, errors in Doppler mea- 
surements are much greater than t he above values and are usually systematic 
( Griffin fc Griffinlll973tlBrownlll990l ). 



2.3.1 Requirements for high-precision Doppler measurements 

Doppler measurements made with spectrometers o r classical spectrographs rarely 

achie ve standard errors better than ~200 m s~^ (e.g. Marcv fc Benit j ll98!^:toua uennov et al 



19911 ). which falls short of what is needed to detect extrasolar planets. To achieve 
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the high Doppler precision necessary for the detection of extrasolar planets, the 
following sources of systematic errors must be addressed: 

1. Motion of the photocenter at the spectrograph slit. In Doppler units, the 
slit width of a spectrograph is typically a few kms~^, meaning that errors 
larger than 1 ms^^ will occur if the photocenter moves away from the slit 
center by a few 10~^ slit widths. In practice, photocenter motions due to 
guiding errors, focus, seeing fluctuations, and atmospheric refraction usually 
amount to '^10^^ slit widths. An efficient solution to this problem is to use 
an optical fiber (often with the addition of a scrambling device) to convey 
and scramble the starlight from the telescope to the spectrograph, thereby 
producing a nearly uniformly illuminated disk at the spectrograph entrance. 

2. Variations of air refractive index and thcrmomechanical flexures. Variations 
of temperature and barometric pressure modify the refractive index of air 
near the grating, causing spurious wavelength shifts similar to Doppler shifts. 
Mechanical instabilities and thermal relaxation also produce non-negligible 
motions of the spectrum relative to the detector, causing additional spu- 
rious shifts and changes in the instrumental profile. For the CORALIE 
spectrograph, a temperature change of 1 K produces a net velocity drift of 
~90 ms~^, while a pressure change of 1 mbar produces a net velocity drift 
of ~300 ms~^. A partial solution to this problem is to stabilize and control 
the entire spectrograph in temperature and pressure. Yet, small wavelength 
shifts cannot be fully avoided and the general way to deal with this problem 
is to use a simultaneous wavelength calibration to monitor - and then correct 
for - these "instrumental drifts" . 

3. Timing of exposure and barycentric correction. Classical radial velocities 
are corrected to the Solar System barycenter. The main contributions to the 
Earth's barycentric motion are the diurnal rotation of the Earth (1-2 ms^^ 
per minute at most) and the Earth's orbital revolution (±30 kms"^ per 
year). Yet, at the ms^"'^ level, the motion of the Sun around the Solar System 
barycenter and the motion of the Earth around the Earth-Moon barycenter 
are also significant. To obtain precise barycentric radial velocities, one thus 
needs precise Solar System ephemeris and one needs to know the photon- 
weighted midpoint of each observation to better than 30 s. 

The above requirements show that the secret to making high-precision Doppler 
measurements is to record the wavelength reference spectrum simultaneously with 
the stellar spectrum and to efficiently scramble the starlight before sending it to 
the spectrograph. In terms of measurements, the observed radial velocity (Kbs) 
must be corrected for the instrumental drift (Mnst, mainly due to points 1 and 2 
above) and for the Earth's barycentric motion (V^Earth, point 3 above) to become 
the stellar barycentric radial velocity T^tar = K>bs — Mnst — ^^Earth used in the 
Keplerian analysis described in Sect. 12.21 Two techniques have successfully been 
developed to accomplish this. If both these techniques use a cross-dispersed echelle 
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Fig. 3. Principle of the two techniques used to obtain high-precision radial velocities. 



spectrograph and track instrumental drifts by means of a simultaneous wavelength 
calibration, they differ fundamentally in their approach and in their design. 

2.3.2 The simultaneous reference technique 

As illustrated on Fig. [3] (lower panel), the simultaneous reference technique involves 
the use of two optical fibers to feed the spectrograph: the so-called object fiber 
which records the starlight, and the so-called reference fiber which records the 
light from a wavelength reference source. In practice, the two fibers are brought 
into the entrance plane of the spectrograph in close proximity to one another, 
separated in the direction perpendicular to the main dispersion. The object and 
the reference spectra are then recorded simultaneously, the two sets of echelle 
orders being distinct and alternate on the detector. T he simultaneous refe rence 
technique was pioneered on the ELODIE spectrograph (iBaranne et al. 1996() , and 



up to now all the spectrographs based on this technique have used as reference 
source a thorium-argon (ThAr) lamp. The method is thus commonly known as 
the simultaneous thorium technique. 

A distinguishing feature of the simultaneous reference technique is that it can 
achieve high radial-velocity precision only if implemented on fiber-fed spectro- 
graphs with high mechanical and thermal stability. The use of optical fibers 
equipped with a scrambling device is indeed essential to reduce instabilities in 
illumination and variations of instrumental profile, both of which cannot be cor- 
rected by the simultaneous calibration. As to the thermal stabilization, it is meant 
to keep the optical paths of the two beams very similar within the spectrograph. 
Under such circumstances, the residual instrumental drifts experienced by the 
two beams are highly correlated and the velocity drift measured by the reference 
channel can be used to correct the object channel. 
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If the simultaneous thorium technique puts considerable demands on the in- 
strumentation, the Doppler analysis is relatively straightforward. Radial velocities 
are traditionally obtained by numerically cross-correlating the o bserved spectra 
with box-shaped, binary (0 and 1 values) templates called masks ( Baranne et al.l 
1996t IPepe et al.ll2002a[ ). Schematically, the cross-correlation function (CCF) is 
constructed by shifting the velocity of the mask by increasing amounts over a 
window roughly centered on the radial velocity of the star (Fig. |4|). The better 
the alignment between the stellar lines and their box-shaped counterparts in the 
mask, the lower the cross-correlation value. The CCF is thus minimal when the 
velocity of the mask perfectly matches the radial velocity of the star (middle panel 
in Fig. 2]). For slowly rotating dwarfs {vsmi < 10 kms~^), the central part of the 
CCF is well approximated by a Gaussian function and the radial velocity of the 
star is measured by the center of the best-fit Gaussian. Note that radial veloci- 
ties obtained through cross-correlation are free from most of the systematic errors 
listed in Sect. 12.3.11 only relative to the stellar mask, which defines the velocity 
zero point. 

The masks used for cross-correlation are of two types: stellar and thorium. 
Stellar spectra are cross-correlated with stellar masks whose nonzero zones corre- 
sponds to the theoretical positions and widths of stellar absorption lines at zero 
velocity. Experience shows that there is no need to use a different mask for each 
star. Common practice is to use a few masks corresponding to the main spec- 
tral subtypes (e.g. G2, K5, and M2). Stellar masks are built either from high- 
resolution, high signal-to-noise observed spectra or from synthetic spectra. Ref- 
erence thorium spectra are cross-c orrelated with a thorium mask built from the 
atlas o f Palmer fc Engleman (1983) and from the updated line list of lLovis &: Pepe 
(2003). 

Stabilized spectrographs searching for planets using the simultaneous tho- 
rium technique include ELODIE (1.93-m Telescope, Haute-Provence, France; dis- 
mounted in 2006), CORALIE (Leonard Euler Telescope, La Silla, Chile), HARPS 
(3.6-m ESQ Telescope, La Silla, Chile), and SOPHIE (1.93-m Telescope, Haute- 
Provence, France). The radial- velocity precision has considerably improved from 
ELODIE to HARPS and is now better than 1 ms"! (Sect. K3\i . 



2.3.3 The gas cell technique 

The basic feature of the gas cell technique is to pass the starlight through a 
cell containing an absorbing medium (the reference source) before entry into the 
spectrograph (Fig. [3l upper panel) . In this way, absorption lines from the refer- 
ence source superimpose on the stellar spectrum, providing a fiducial wavelength 
scale that experiences the same instrumental shifts and distortions as the stel- 
lar spectrum. Since the gas cell is at rest relative to the observatory, spurious 
instrumental shifts are measured by the wavelength shift of the reference lines. 
The reference spectrum also provides a specification of the instrumental profile 
(IP) at each position on the detector, allowing for the measurement and correc- 
tion of IP variations. The first application of the gas cell technique to planet 
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Fig. 4. Illustration of the construction of the cross-correlation function. Diagrams on 
the left represent the stellar spectrum (dashed lines) and the binary mask (solid lines, 
transmission zones depicted as hatched areas). Diagrams on the right show the result of 
the cross-correlation process. Courtesy of Claudio Melo. 



sear ches was made bv ICampbell et al.l (jl988l ). who used a hydrogen fluoride (HF) 
cell (ICampbell fc Walked I igTolT At present, all planet search programs based on 
this technique use an iodine cell and the technique is commonly referred to as the 
iodine cell technique. 

Instrumentally, the iodine cell technique is easily implemented on any existing 
slit spectrograph. The main complication of the method resides in the Doppler 
analysis. In practice, spectra taken through the iodine cell are broken up into 
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Wavelength (A) 

Fig. 5. Illustration of the modeling process for the iodine cell technique. Top: transmis- 
sion function of the iodine cell (Tja) as measured by the FTS. Second: a template stellar 
spectrum (proxy for Is). Third: the composite spectrum of the same star, as observed 
through the cell (dots) and as modeled (solid line). Bottom: ten times the difference 
between the model and the observations. Figure from iButler et al.l il99dl ). 

several hundred "chunks" of length ^^2 A. On each chunk the composite spectrum 
is modeled as 

lobs (A) = k [Ti, (A)/, (A + AA)] * IP (2.8) 

where Is is the intrinsic stellar spectrum, T/^ is the transmission function of the 
iodine cell, IP is the in situ instantaneous instrumental profile, and the constant 
fc is a normalization factor. The wavelength shift, AA, is the topocentric Doppler 
shift of the star. The barycentric radial velocity is obtained by correcting this 
topocentric Doppler shift for the Earth's barycentric motion and by converting the 
result into a radial velocity using an elaborated version of the Doppler formula. 
Again, the barycentric radial velocity obtained in this way is not an absolute 
velocity; it is precise only relative to a stellar template, which defines an arbitrary 
velocity zero point. 

As shown by Eg. 12.81 the modeling process requires two input functions, and 
Tj^, plus the knowledge of the instrumental profile. The transmission function of 
the cell is measured directly by obtaining a spectrum of the cell with a Fourier- 
Transform Spectrometer (FTS). The intrinsic stellar spectrum is more difficult to 
obtain since it cannot be measured directly. is usually generated by taking a 
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high signal-to-noise spectrum of the program star without the iodine cell in place 
(which gives Ig * IP) and by deconvolving this spectrum from the instrumental 
profile. The instrumental profile cannot be measured directly either. It is com- 
monly generated by comparing observations of a hot, rapidly rotating star (which 
basically gives Tj^ * IP) with the reference FTS iodine spectrum. 

The whole modeling process is illustrated on Fig. [5] and typically requires 13 
parameters: 2 for the wavelength scale, 1 for the topocentric Doppler shift of the 
star, and 10 for the description of the instrumental profile. These parameters are 
determined by comparing a series of mod els to the observ e d spe c trum through 
a leas t-squa res fit. We r e fer th e reader to Marcv fc Butler (1993), Valenti et al 



( 1995h . and iButler eTaP ('1996') 'or complementary information on the practical 
aspects of the iodine cell technique. 

Echelle spectrographs equipped with an iodine cell and actively used for planet 
searches include HAMILTON (Shane/CAT, Lick, USA), HIRES (KECK, Mauna 
Kea, USA), UCLES (AAT, Siding Spring, Australia), HRS (HET, Davis Moun- 
tains, USA), and UVES (VLT, Paranal, Chile). The iodine cell technique has 
long demonstrated a precision of 3 ms" (iButler et al.lll996h . but the precision 
seems to have not much improved since then. Residuals arou nd the best plane- 
tary orbits published so far indicate a precision of ^2.5 ms~^ (jButler et al. I I2OO6I: 
Fischer et aLll2006t Imight et al.ll2008f) . 



2.4 Limitations on Doppler spectroscopy 

Limitations on Doppler measurements can be classified into three broad categories: 
photon count, technical, and astrophysical. We describe here photon count and as- 
trophysical limitations. Technical limitations related to the simultaneous thorium 
technique will be discussed in Sect. 14.3] 



2.4.1 Photon comit limitations 



The ultimate limit to the attainable Dopple r prec ision is that set by photon- 
counting statistics. As proposed by IConnej ( 1985i ). the uncertainty related to 



photon noise can be quantified by a quality factor, Q, which is a sole function 
of the line profile in the case of pure photon noise. For spectral types between 
F2V and K7V, Q incre ases towards the blue and Q is the highest for a K5V 



star (jBouchv et al.l 120011 ). These two trends reflect the evolution in the strength 
and width of metal lines along the spectral sequence (late-type stars have more 
"peaky" lines) and the fact that these lines are more numerous between 4000 
and 4500 A than between 6000 and 6800 A. Q is also sensitive to rotational line 
broadening, decreasing as 1/v sini for v sini > 6 kms~^. For these reasons, early- 
type stars (^F6) with few and broad absorption lines do not permit high-precision 
Doppler measurements, whereas slowly rotating K dwarfs are ideally suited for 
high-precision Doppler programs. 
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Fig. 6. Illustration of the effect of a starspot on the profile of the cross-correlation 
function and on the determination of its centroid (radial velocity). The effect has been 
exaggerated for illustrating purposes. 



2.4.2 Astrophysical limitations 



Stellar activity Magnetic phenomena at the surface of solar-type stars induce 
radial- velocity variations through the t emporal and spatial evolution of spots , 



plages, and convective inhomogeneities ( Saar &: Donahue! 1997 : Saar et al. 1998 ) 



Figure [6] illustrates how the rotation of a starspot modifies the line profiles - and 
equivalently the CCF profile - affecting the determination of the radial velocity. 
When the starspot pattern is long-lived, variations in line asymmetry are modu- 
lated by the rotational period of the s tar and can mimic a planetary signal (e.g. 
Queloz et al.l 12001 : Bonfils et al. 2007 ). When the star is observed longer than 
the typical lifetime of starspots, the signal becomes incoherent and is detected as 
radial-velocity "noise" . Although less confusing, this is still annoying since it may 
inhibit the detection of planetary signals of lower amplitude. 

In practice, radial- velocity perturbations related to stellar activity are referred 
to as "stellar jitter" . Stellar jitter de pends on effective temperature, stellar activity, 



and projected rotational velocity (jSaar et al.l 119981 : iSantos et al.l l2000bt iWright 



200i), but these dependencies cannot be modeled in detail yet. Typical values of 
stellar jitter are < 5 ms~^ for slow ly rotating, chromospherically quiet G-K dwarfs 
Santos et al.ll2000b| :lwigh^l2005l) and < 50 ms'^ for F5-M2 young, active dwarfs 



Paulson et al.ll2004D . The above values explain why high- precision Doppler planet 



searches select their targets among old, slowly rotating inactive stars. 

To quantify the activity level of their (potential) targets, Doppler planet searches 
traditionally use the parameter i?HK > which represents the fraction of a star's bolo- 



metric flux emitted by the chromosphere in the Ca ii H and K lines (jNoves et al 



1984 . Since this chromospheric emission is closely related to the surface magnetic 
flux, a high R'^^ value is an indication that a star may exhibit significant activity- 
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Fig. 7. Sequences of HARPS radial-velocity measurements on four solar-type stars. 
Solar-like oscillations are clearly visible on each of these sequences. 



related velocity variations. But a high R'^^ index does not prove that a given 
Doppler signal is of stellar origin. To prove this, one can search for a periodic 
modulation in the Ca ii flux, i n the bisector of the cross-correlation function, or in 
photo metric observations (e.g. Queloz et al. 200ll : lDesort et al1l2007 : Bonfils et al.l 



20071 ) . These diagnostics are pretty effective at identifying activity-related velocity 
variations in G, K, and M dwarfs. 



Solar-like oscillations The Sun oscillates in many low-amplitude modes simul- 
taneously, each mode being characterized by a maximum amplitude ^20 cms~^ 
and a period between 3 and 15 minutes. This so-called five-minute oscillation is ex- 
cited stochastically by turbulent convective motions near the solar surface. Main- 
sequence stars with significant outer convection zones should all exhibit similar 
pulsations, called solar-like oscillations in this more general context. Recent ob- 
servations have shown that this is indeed the case (Fig. [7]). It should be emphasized 
that while individual modes have typical amplitudes in the range 10-50 cms~^ for 
main-sequence stars, their interference can produce radial-velocity variations up 
to about 10 times these values on time scales of minutes. Although these oscilla- 
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tions constitute a very interesting tool to probe the structure of stellar interiors 
(this is asteroseismology; see the contribution by Vauclair), they are essentially 
"an annoying source of noise" when viewed from the side of very high precision 
planet searches. As explained in Sect. 14.31 an adequate observing strategy allows 
fortunately to largely average out this noise, reducing its impact on the search for 
low-mass planets. 



Contamination effects Slit spectrographs used for planet searches have typi- 
cal on-sky slit widths of 0.5-1", while fiber spectrographs have typical on-sky fiber 
diameters of 1-3". The light from an object close to the scientific target may thus 
also fall within the slit/fiber, contaminating the target's spectrum and altering 
the determination of the target's velocity. For this reason, known binary stars 
closer than 2-6" are genera lly rejected from Doppler planet search programs (see 
Eggenberger fc Udrv 20071 for details). Yet, some close star systems are present 



in the samples of large-scale Doppler surveys and it may happen that a star turns 
out to be a binary, or that a presumed binary system turns out to be triple. The 
latter case is of particular concern because unrecognized triple systems may under 
some circumstances i nimic the radial- ve l ocity signature of a planetary companion 
to a single star fe.g. ISantos et al.ll2002l : IZucker et al.l 120031 ). Unrec ognized triple 
systems can be ident ified indirectly through bisector analyses (e.g. ^Sant os et al 



2002.; Eggenberger et a l. 200 3) and direct l y usin g two-dimensional correlation (e.g. 
Zucker et al.l l2003r Eggenberger fc Udrv 2007f ). Due to the strong bias against 



(moderately) close binaries in the samples of regular Doppler planet searches, con- 
tamination effects remain marginal and do not significantly affect planet searches. 



3 Statistical properties of planetary systems with giant planets 

The number of planetary companions detected through Doppler spectroscopy has 
increased exponentially since 1995, yielding an abundant observational material to 
characterize the properties of giant planets orbiting within 3 AU from their host 
stars. These data have been used extensively to study the distributions of planet 
properties and to q uantify the iiicidence of giant planets among different popu - 
lations of stars (e.g. iMarcv et al1l2005at lUdrv et allboOTt lUdrv fc San"to^l2007f ). 
As these statistical features are thought to bear fossil imprints of planet forma- 
tion and evolution processes, they provide a crucial link between observation and 
theory. This link has recently been strengthened by the new ability of core accre- 
tion models to generate synthetic planet po pulations, which can be compared with 
the observed sample of extrasolar planets ( Ida fc Lin 2004at lAlibert et al. 2005t 
Mordasini et al. 2009, in prep.). 



3.1 Doppler surveys and their samples 

Radial- velocity planet search programs have traditionally selected their targets to 
optimize the achievable Doppler precision, favoring bright, chromospherically in- 
active main-sequence stars with spectral types between ~F8 and ~M0. Among 
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these "classical" programs, those which have most significantly contri buted to the 
prese nt discoveries are: the ELODIE Planet Search (~330 G-K dwarfs :|Perrier et al. 
20031 ). the California and Carnegie Planet Search ('^llOO F- G-K-M dwarfs; Marcv et al. 



2005bH . the CORALIE Planet Search (~1650 G-K dwar fs: lUdrv et al[l2000h . the 
Anglo Austrahan Planet Search (--230 F-G-K-M stars: iTinney et al.ll2005l ). and 
the HARPS Planet Search (-1500 G-K-M dwarfs: iMavor et al.feoo'i r To this list 
we can add the so- called "metallicity-biased" surveys, N2K (—2000 F-G-K dwarfs; 
Fischer et all2005h and ELODIli (-1200 G-K dwarfs: Ida Silva et al.ll2006h . which 
aim at detecting short-period giant planets orbiting metal-rich stars. Omitting the 
metallicity-biased surveys, these programs collectively monitor about 3000 F-G-K 
dwarfs of the solar neighborhood. 

Planet searc hes targeting M dwarfs were originally lim ited by the faintness 
of their targets (jPelfosse et al. I ll998l: iMarcv &: Butled Il998h . Recent instrumen- 
tal developments coupled with the use of larger telescopes have considerably im- 
proved the situation, all owing the inclusion of a significant number of M dwarfs 
in the present surv eys (|Endl et all 120031 : iButler et al] [20041 : [Bonfils et all 120051 : 
Bouchy et al.l 120071 ). In total, about 300 nearby M dwarfs are now being moni- 



tored. 

As explained in Sect. [2^ main-sequence stars earlier than — F6 are not suited 
for high-precision Doppler measurements and classical surveys avoid them accord- 
ingly. A recent attempt has been made to extr act precise radial ve locities for 



A and F dwarfs following a dedicated approach (jCalland et al.l 120051 ) . but even 



with this method planet searches remain limited by stellar activity. An alterna- 
tive approach to study planet occurrence among stars more massive than the Sun 
consists in observing these stars when they have evolved off the main-sequence, be- 
coming G-K giants or subgiants. Doppler planet searches around G-K (sub)giants 
have gradually expande d during the past five years and some hundred stars are 
now regularly observed |Frink et al ]l2002l : ISetiawan et al.ll2003l: iHatzes et al.ll2005t 
Johnson et alll2006HLovis fc Mavoil2007l iNiedzielski et al.ll2007tlSato et al.ll2005l ). 



A notable built-in bias affecting classical Doppler planet searches is the avoid- 
ance of binaries closer than 2-6" fSect. Due to this discrimination, current 
data only provide sparse information on the suitability of J;200 AU binaries for 
planetary systems. But again the situation is improving, and radial-velocity planet 
searches have recently been extended to spectroscopic and modera tely close visual 
binaries (e.g. lEggenberger fc Udrv 2007 : Muterspaugh et al1l2007t ). 

Turning to statistical studies, it should be noted that rigorous analyses includ- 
ing the calculation of detectio n thresholds and surv e y completeness are available 
only for the oldest prog rams (Gumming et al. '1999I: fandl et al.ll2002l : iNaef et al 



2005t Wittenmyer et al. 2006 : Gummin g et al, 20081) , and with the exception of 



^BB^ML-^^^^^^^^^^S^r^ ^^^^^^m^ I - 1; - ■ ^-T^^^^^ . ^ . r^~. ^ — — „— .-^ TV . jM/ ; ■ JT 

Gumming et al.l ( 20081 ) these analyses contain too few planet detections to study 
the distributions of planet properties. Therefore, even though they have not been 
rigorously corrected for selection effects yet, the two surveys best suited for sta- 



2 Now replaced by a SOPHIE program llBouchv et al.ll2007ri . 
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tistical analyses are the CORALIE survey with its HARPS extension and the 
combined Lick+Keck+AAT survey. 



3.2 Occurrence rate of giant planets 

At present, true occurrence rates are available for the ELODIE and for the Keck 
programs. Results from the ELODIE survey indicate that the fraction of stars 
with a planet more massive than 0.5 Mj^p is 0.7 ± 0.5 % for P < 5 days (a < 0.06 
AU) and 7.3 ± 1.5% for P < 3900 days (a < 4.8 AU) (|Naef et al.ll2005l ). For the 
same mass and period ranges, the Keck survey obtains 0.65±0.40% and 8.6±1.3% 
(jCumming et al.ll2008l ). 

It should be emphasized that according to planet formation models, giant plan- 
ets reside preferentially beyond ~3 AU, whil e the overall planet population is 
domi nated (in number) by low-mass planets ( Ida fc LinI 2004a : Mordasini et al 



20071 ). Taking into account the selection effects inherent to Doppler spectroscopy, 
the ELODIE and Keck programs coul d detect only a small f raction of the overall 
planet population (^6% according to iMordasini et al.l 120071 ). The above figures 
thus lend support to the hypothesis that planets are common around solar-type 
stars. 



3.3 Mass, period, and eccentricity distributions 
3.3.1 Mass distribution 



As explained in Sect. 12.21 Doppler spectroscopy does not yield the true mass of 
planetary companions but only the pro duct rriri sini, which is a lower limit on the 
planet mass. Nonetheless, as shown bv I Jorissen et al. ( 200ll ) . the distribution of 
minimum masses is a good proxy for the mass distributioio 

The distribution of minimum masses is shown in Fig. [H While most of the 
detected planets have minimum masses below 5 Mjup, the distribution exhibits a 
long tail towards minimum masses larger than 10 Mjup. For orbital periods less 
than a decade, there is a s carcity of companions w ith (minimum) masses between 
~10 and ~100 Mj^p (e.g. IZucker fc Mazehl 120011) . This gap in the mass distri- 
bution - called the brown dwarf desert - constitutes the most obvious evidence 
that extrasolar planets and stellar companions belong to two distinct populations 
which formed through two different mechanisms. The two populations are not 
completely disconnected, however, as some overlap seems to exist between ~10 
and ~20 Mjup. Distinguishing massive planets from low- mass brown dwarfs on 
the sole basis of their (minimum) mass may thus be hazardous, preventing at 
present a determination of the upper limit to the mass of a planet. 

The low-mass side of the minimum mass distribution is strongly affected by 
selection effects below the mass of Saturn (0.3 Mjup), but despite this incomplete- 



^For the sake of conciseness, we often use the terms "mass" and "mass distribution" instead 
of "minimum mass" and "minimum mass distribution". However, we make the difference or use 
the additional term "true mass" when the difference is important. 
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MjSini [Mj„p] log( M^sini ) [Mj^^] 



Fig. 8. Distribution of planet minimum masses in linear {left) and logarithmic (right) 
scales. The double-hatched histogram represents the planets detected with CORALIE 
(left) or with HARPS (right). 



ness the distribution seems to turn upwards near 0.1 Mjup. This low-mass part of 
the distribution will be discussed in Sect. 14.2."^ 

3.3.2 Period distribution 

The period distribution is shown in Fig. [9l At short periods, there is a pile-up of 
planets with periods between 2 and 5 days. These are the so-called hot Jupiters like 
51 Pegb. The detection of Jovian planets orbiting their stars in just a few days has 
been a major surprise since prior to this discovery theory predicted that giant plan- 
ets would form - and by implication exist - only beyond a fe w AU. Although not 



comp letely excluded, forming these planets in situ is difficult (jBodenheimer et al 



200011 and the first detections of hot Jupiters have posed a serious challenge to 



planet formation theories. According to the revised models, short-period giant 
planets likely formed beyond a few AU and the n migrated inwards as a result of 
their tidal interaction with the gaseous disk fe.g. lLin et ai1ll996l : llda fc Linll2004at 



LlVIordasini ct al. 2007: lArmitagel 2007 1. The discovery of close-in giant planets like 



51 Pegb has thus revived the notion of planetary migration (see the contributions 
by Terquem and Ida) and there is now little doubt that many extrasolar planets 
have formed at locations which do not correspond to those where they are observed 
today. 

The origin and even the existence of the pile-up of planetary companions with 
periods close to 3 days are still being debated. If hot Jupiters have undergone 
extensive migration within a gaseous disk, there is a priori little reason why they 
should have stopped at small radii rather than plunge into their star. Some 
mechanisms have thus been proposed to stop orbital migration and to park gi- 
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log(Period) [days] 

Fig. 9. Period distribution for extrasolar planets. The hatched histogram represents light 
planets with m^sini < 0.75 Mjup, while the filled histogram represents Neptune- mass 
planets with mp sin i < 21 M^. 



ant planets on close orbits: the exist ence of a central cavity in the disk (e.g. 
Lin et al.[ 1996t Kuchner fc Lecaj 2002), tidal intera ctions with the spinning host 



star (e.g. [ Tri lling et al. 1998; Patz old fc Raueiil2002n . or Roche lobe overflow (e.g. 
Trilling et al.. 1998: Gu et al. 2003) . Alternative scenarios for the fo rmation of h ot 



Jupit ers involve resonant interactions with a disk of planetes imals (|Murray et al 
19981). or tidal ci rcularization of highly eccentric orbits (e.g. Ford fc Rasiol 2007t 



Fabrvckv fc Tremaine 2007) . 



At longer periods, the number of planets remains fairly constant between ^^7 
and ~250 days, and then increases abruptly. The time baseline of ongoing pro- 
grams is still too short to tell whether the period distribution (in logarithmic 
scale) remains flat from ^300 days to ~4 years, or increases continuously beyond 
~250 days. In either case, current data indicate that there exists a large reservoir 
of giant planets with periods above 5 years (a ^ 3 AU) which we are just beginning 
to detect. This shows that the main limitation at the moment for the detection of 
analogs of our Solar System is not the Doppler precision but the duration of the 
surveys. 



3.3.3 Period-eccentricity diagram 

Figure [10] shows the period-eccentricity diagram for both extrasolar planets and 
stellar binaries with solar- type primaries. This plot makes it clear that most 
of the extrasolar planets have anomalous eccentricities compared with those of 
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Fig. 10. Period-eccentricity diagram for extrasolar planets (open pentagons) and for 
stellar companions to solar-type stars (black dots) . The Earth and the four giant planets 
of the Solar System are also indicated (Earth and starry symbols at the bottom). 



Jupiter and Saturn (0.048 and 0.056, respectively). Planets with pe riods below 
^5 d ays tend to have small eccentricities due to tidal dissipation (e.g. iRasio et al.l 
1996f) . Tidal circularization becomes ineffective at longer periods and planetary 



eccentricities then span almost all the allowable range with a median value of 0.29. 

These high planetary eccentricities are very puzzling since formation in a disk is 
thought to yield almost circular orbits. If it is so indeed, some additional processes 
must drive the planets out of their initially circular orbits . Numerous mechan isms 
have been proposed to excite planetary eccentricities (see iFord fc Rasidl2007l for a 
review and references), including strong gravitational scattering between planets, 
secular interactions with a distant companion in a highly inclined orbit, or interac- 
tions of orbital migration with mean motion resonances (see also Sect. 13. 4| ). While 
these three mechanisms are all supported by some observational evidence, none 
can reproduce the overall eccentricity distribution. It is thus likely that several 
mechanisms contribute in determining the eccentricity distribution. 

Contrary to the mass distribution, the eccentricity distribution of extrasolar 
planets is not radically different from th at of stellar companio ns to solar-type 



stars, though small differences exist (e.g. Halbwachs et al. 20051 ). This observa- 



tion has sometimes been considered as an indication that extrasolar planets and 
stellar companions belong to the same population and share a common origin 
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( Stepinski fc Black 20001 ). Actually, it is more probable that the formation mecha- 
nisms were distinct, but that the orbital eccentricities were excited through similar 
mechanisms. 



3.4 Multiple planet systems 

Among the 253 planet-host stars with Doppler measurements, 30 are orbited by 
multiple planets, th e most prolific systems bein g presently 55 Cn c with 5 planets 
([Fischer et al.ll2008[) and fiAra with 4 planets (jPepe et ahlbOOTl ). The observed 
fraction of multiple planet systems is 12%, meaning that the probability of finding a 
second planet is enhanced by a factor of about two with respect to the probability of 
finding the first planet. This multiplicity rate is certainly a lower limit as multiple 
planet systems require more Doppler observations spanning a longer time baseline 
to extract the different Keplerian signals. In the surveys with the longest duration 
the fraction of multiple systems ranges f rom 25% to 50% ( Fischer fc Valentil 20051 : 
Udrv fc Santo^l2007t IWright et al.ll2007l ). These results suggest that like the Sun, 
most solar-type stars are orbited by systems of planets. 

The presence of two (or several) planets interacting one with another is in- 
teresting because it increases our potential ability to characterize the system's 
orbital configuration and formation history. For instance, strong mutual interac- 
tions between two planets a llow to (partially) remove the sin i dege neracy, yielding 



true planetary masses (e.g. IChiang et al]l200ll : iRivera et al1l2005f ). From a more 



specific point of view, the existence of pairs of planets locked in low-order mean 
motion resonances (i.e. with commensurable periods such as GJ 876 b a nd c) indi- 



cates tha t some planets experienced smooth convergent migration (e.g. |Leg|2004; 
Kley et al. 2005). Alternatively, the current orbital configuration of the three gi- 



ant planets orbiting Ups And suggests that planet-planet scattering occurred in 
this system (|Ford et al.ll2005l ). 



3.5 Properties of the host stars 

Planet formation being a by-product of star formation, some stellar characteristics 
may provide important information about the conditions leading to planet forma- 
tion. Recent work has shown that the occurrence rate and the properties of giant 
planets depend on the mass, metallicity, and multiplicity status of the parent star. 



3.5.1 Stellar metallicity 

Not long after the discovery of the first extrasolar planets, it was noticed that 
planet-bearing stars were systemat ically metal-rich compared to the Sun (e.g. 



Gonzalez 1997 : Santos et al. l2000ah . Recent studies have shown that planet oc- 
currence is a rising function of the parent star metallicity (Fig. Ilip . Specifically, 
25% of the stars with twice the metal content of the Sun ([ Fe/H]= 0.3) harbor 



a giant planet, against only ~3% for solar-metallicity stars ( Santos et al.ll2004bl : 



Fischer fc Valentil 120051 ). This observation has motivated the metallicity-biased 



surveys mentioned in Sect. 13.11 
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Fig. 11. Percentage of planet- host stars as a function of stellar metallicity. The dashed 
line shows the results from lSantos et all (|2004lJ ) and the solid line shows the results from 
iFischer fc Valentil l|2005h . 



Three explanations have been proposed to understand this correlation: (1) an 
enhanced metallicity in the protoplanetary disk could be more conductive to planet 
formation, (2) the infall of metal-rich planetary material onto the host stars could 
have enriched their outer layers, or (3) there might be an "orbital period bia s" due 
to the dependence of migration rates on metallicity ( Livio fc Pringle 2003 h or to 
the i nward shift of the o ptimum region for giant planet formation in metal-poor 
stars ( Pinotti et al.ll2005 ). Some observat ional evidence for the last two hypot heses 
have been reported in the literature (e.g. Ilsraelian et aL 2001 : Sozzetti 2004 ). but 
most analyses suggest that the overall metallicity excess observed in planet-host 
stars has a pr imordial origin and reflects the high metallicity in the protoplane- 



tary disk (e.g.lPinsonneault et al.l l200ll: iLivio fc Pri^^l2003l : ISantos et al.H2003t 



Fischer fc Valentil l2005f ) . " This conclusion is c orroborated by core accretion mod 



els, which reproduce the observed correlation ( Ida fc Linll2004lj: Benz et al. 2006f ) 



Disk instability being apparently insensitive to metallicity (|Bossll2002l ) , the metal 



rich nature of planet-host stars is commonly considered as a major evidence that 
most of the giant planets we observe today formed through core accretion. 



3.5.2 Stellar mass 



The occurrence rate of Jovian planets appe ars to correlate w i th stellar mass, gi- 
ant planets being rarer around M dwarfs et all 120061: iButler et al.l 120061: 



Bonfils et al.ll2007| ) and more frequent around intermediate-mass stars (|Lovis fc Mavor 
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20071: 1 Johnson et al.ll2007f ). Another quantity which corr elates with stellar mas s 
is the "average mass of planetary systems" as defined by iLovis fc Mayor ( 2007 ). 
This second trend may simply be a consequence of the first correlation, but it could 
also indicate that planetary systems tend to be more massive around more mas- 
sive stars. Under the assumption that the mass of protoplanetary disks scales with 
stellar mass - which r emains to be confirmed - both trends a.re supported by the 
core accretion model (Laughlin et al. 2004; Ida fc Lin 20051 : Kennedy fc Kenvon 
20081: iBenz et~al]l2008h . 

Remarkably, G-K (sub)giants hosting gian t planets seem to lack the strong 
meta l-rich tendency seen in F-G-K dwarfs (e.g. iPasauini et aL 2007; Takeda et al 
2008f) . It is thus likely that both metallicity and stellar mass impact the occur- 



rence and the properties of giant planets. Larger samples of planets orbiting low- 
and intermediate-mass stars will be needed to confirm this observation and to 
disentangle the two effects. 



3.5.3 Stellar multiplicity 

Contrary to some expectations and despite the strong bias against their detec- 
tion, giant planets have been found in various types of binaries and triple sys- 
tems, including spectroscopic binaries with projected separations of ~20 A U (e.g. 
Eggenberger et al.ll2004 iRaghavan et aLluOOa IPesidera &: Barbieril 120071 ). Inter- 
estingly, some of these planets seem to possess distinctive characteristics when 
compared to planets orbiti ng single stars: (1) the most massive short-period plan- 
ets ar e all found in binaries ( Udry et al. Il2002l : IZucker fc Mazehll2002l : lEggenberger et al 
2004 ) ■ (2) the minimum period for a significant eccentricity seems larger for 
planets in binaries (■-^4 days) than for planets around single stars (~5 days; 
Eggenberger et al. 2004[) . and (3) the four plane ts with an eccentric ity larger than 
0.8 all have a stellar or brown dwarf companion ( Tamuz et al. 20081) . These trends 



may indicate that migration and postformation evolution proceed differently in 
some ty pes of bina r ies than around single stars, as suggested by som e theoretical 
studies (|Kley||2000l: IWu fc Murray||2003HFabrvcky fc Tremainel[2007l) . 



4 The road to Earth twins 



One of the main drivers of extrasolar planet searches is to understand the origin 
and the degree of uniqueness of a life-bearing planet like the Earth. There are thus 
strong motivations to search for rocky planets orbiting within the habitable zones 
of their stars, that is within the circ umstellar reg ion where a terrestrial planet can 



hold Uquid water on its surface (e.g. lKasting et al.„ 19931. Accord ing to theoretical 
mode ls, such planets should be numerous re.g. Ilda fc Linll2004al : iMordasini et al.l 
2007 ) . The main obstacle to their discovery through Doppler spectroscopy is their 



20071) . The main obstacle to their discovery through Doppler spectroscopy is their 
low mass, which translates into an extremely low-amplitude radial-velocity signal 
(9 cms~^ for the Earth). 

The detection of low-mass planets being closely related to the ultimate precision 
achieved by Doppler surveys, the first Neptune- mass planets (10 M0 < mpSini < 
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25 Mffl) were found only four years ago (jButler et al.ll2004 iMcArthur at al.|[2004 : 



Santos at al. 2004a[) . The first observation of a so-called super-Ear th (2 M 



< 

mpSini < 10 M^) was announced one year later (jRivera et al.ll2005D and opened 
the way to a series of discoveries. This new step forward in planet searches was 
made possible first by the development of a new generation of very precise instru- 
ments whose prototype is the HARPS spectrograph, and second by the application 
of a careful observing strategy to improve detection capabilities and minimize the 
impact of stellar noise. 



4.1 HARPS, a new-generation instrument to search for low-mass planets 

In 1998 the European Southern Observatory (ESO) issued an announcement of 
opportunity for the procurement of an instrument dedicated to the search for ex- 
trasolar planet at the unequaled precision of 1 ms~^. This led to the development 
of the HARPS spectrograph, which was manufactured by the HARPS Consor- 
tiunfl. HARPS is mounted on the ESO 3.6-m telescope at La Silla Observatory 
(Chile) and it has been available to the astronomical community since October 
2003 

HARPS (High Accuracy Radial velocity Planet Searcher) is a fiber-fed, cross- 
dispersed echelle spectrogr aph covering the spectral range from 380 to 690 nm with 



a resolution R = 115,000 (jPepe et all 12002b: iMavor et al 
2004^. Like its two predecessors, ELODIE (jBaranne et al 



20031: iRupprecht et al 



19961 ) and CORALIE 



( Queloz et al. I I2OOOI) . lARPS uses the simultaneous thorium technique. But un- 
like its two predecessors, HARPS is placed inside a vacuum vessel, the operating 
pressure being kept below 0.01 mbar. The vacuum vessel is in turn located in a 
thermally stabilized enclosure and this two-stage insulation ensures a short-term 
(1 night) temperature stability of a few millidegrees. Thanks to this high environ- 
mental stability, instrumental drifts never exceed 1 ms^^ during the night, making 
HARPS the most stable spectrograph in the world. 

In exchange for the delivery of the instrument to ESO, the HARPS Consortium 
received Guaranteed Time Observations (GTO; 500 observing nights over 5 years). 
Since July 2003, the co nsortium has used this time to carry out a comprehensive 
planet search program ( Mavor et al. 2003f ). More than 60% of the GTO time is 
devoted to the search for low-mass planets around '^520 G-K-M dwarfs. 



4.2 Observational results on low-mass planets 

The present sample of low- mass planets (defined here as planets with Wp sin i < 



25 Mq) comprises 18 objects (|Mavor fc Udrvll200 8' and references therein: iForveille et al ]l2008l) 
We present below some noteworthy detections made with HARPS in the context 



*The HARPS Consortium is composed of the Geneva Observatory (leading institute), the 
Observatoire de Haute-Provence, ESO, the Physikahsches Institut der Universitat Bern, and the 
Service d'Aeronomie du CNRS. 
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Fig. 12. Radial velocity measurements of /i Arae obtained during the asteroseismology 
campaign. White circles correspond to night averages. The dashed line shows the Kep- 
lerian orbit of the Neptune-mass planet confirmed by additional Doppler measurements. 
Figure from iBouchv et all (|2005al ). 



of the GTO program. We then discuss a few interesting statistical trends that are 
emerging among this new population of low-mass planets. 

4.2.1 Noteworthy HARPS discoveries 

A Neptune-mass planet around ^ Arae The lO.S-M^ planet around /i Arae 
(/i Araec with a pe riod of 9.6 da ys) was found serendi pitously during an asteroseis- 



mology campaign (ISantos et al . 2004a: iBouchv et al., 2005a'). Prior t o this detec 



tion, /X Arae was known to harbor a giant planet in a 743-day orbit (IButler et al 
2001 ) and a possible additional planetary companion in a wider orbit (|Tones et al 



20021 ). The continuous monitoring of /i Arae within th e HARPS GTO program 
led to the discovery of a fourth planet (jPepe et al.l2007t) . making this star host to 
a rich and complex planetary system. 

As shown on Fig. [121 the amplitude of the radial- velocity signal of ^ Arae c is 
similar to the peak-to-peak radial-velocity variation induced by solar-like oscilla- 
tions. The detection of this tiny planetary signal thus strongly benefitted from 
the high measurement density of the asteroseismology campaign. Nonetheless, 
follow-up observations have shown that adopting an adequate observational strat- 
egy fSect. l4?3)) . planets such as /i Araec are within the detection capabilities of the 
HARPS GTO program. 

The quasi-simultaneous discovery of Neptune-mass planets around 55 Cue, 
GJ 436 and Arae raised the question of their origin and constitution. Are they 
sub- or near-critical planetary cores that failed to accrete large amounts of nebular 
gas (see Sect.[T]), or are they the remnant of gaseous giant planets which lost a la rge 



fraction of their atmosphere through evaporation? While ISantos et ahl (|2004a[ ) fa- 



vor an essentially rocky composition for /x Araec, the gas giant remnant and the 
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ice- ro ck composition cannot be excluded (jBaraffe et al.l 12006 
20051) . 



Brunini fc Cioncol 



A trio of Neptunes around HD 69830 Two years ago, iLovis et al. I (l2006l) 

reported the discovery of a planetary system made of three Neptune-mass objects 
(projected masses of 10.2, 11.8, and 18.1 M®; periods of 8.7, 31.6, and 197 days) 
around HD 69830 (Fig. [T31). This system possess two interesting characteristics: 
it is composed exclusively o f low-mass planets, and it seems to har bor a massive 
astero id belt within 1 AU ( Beichman et al.ll2005f ). According to Alibert et al 



( 2006f ). the iimermost planet is essentially a rocky core, while the two outermost 
planets are rocky cores surrounded by a shell of fluid water and a gaseous envelope. 

From the technical perspective, the values of the radial- velocity semiamplitudes 
(2.2, 2.7 and 3.5 ms^^) make the three planetary signals undetectable with most 
of today's spectrographs. Considering only the most recent measurements made 
with an improved observing strategy, the root-mean-square (rms) of the residuals 
around the best-fit solution is 0.64 ms~^ and this value decreases to 0.35 ms^^ 
(bottom right panel of Fig. [T3|) when the observations are binned to yield one 
measurement per observing run (typical run duration of 10 days). 



Two habitable super-Earths around G1581? The M4 dwarf G1581 also 
harbors a very interesting planetary system made of a hot Neptune {nip sin i = 
16 Mq, period of 5.4 davs: Bonfils et al.,2005.) and two sup er-Earths {mpS'mi = 



5 and 8 M0, periods of 13 and 83 days; lUdrv et al.l I2007D . These two super 



Earths lie close to the inner and outer edges of the habitable zone of G1581, 
raising the question of their habitability. Models of planetary at mospheres indicate 
that planet c is prob ably too close to G1581 to be habitable (jSelsis et al.l 12007 



von Bloh et al. 20071) . Habitability prospects are much better for planet d, which 



may be an interesting target for future space missions aimed at searching for life 
signatures in the atmospheres of low-mass planets. 

This system illustrates two advantages M dwarfs possess over solar-type stars: 
(1) the smaller mass of M dwarfs makes the detection of low-mass planets eas- 
ier fEq. 12. 5[) . and (2) th e habitable zone oc curs farther in around M dwarfs than 
around solar- type stars (j Kasting et al. 1993). These differences imply that detect- 
ing Earth-mass planets in the habitable zones of M dwarfs is much less demanding 
in terms of Doppler precision than detecting the same planets around solar-type 
stars. The planetary system orbiting Gl 581 shows that detecting potentially habit- 
able planets around M dwarfs is within the reach of current high-precision Doppler 
surveys. 



A trio of super-Earths around HD 40307 The latest highlight from the 
HARPS G TO program is th e planetary system discovered around the K dwarf 



HD 40307 ( Mavor et al. 20081 ). This system comprises three super-Earths (mini- 



mum masses of 4.2, 6.8 and 9.2 M^; periods of 4.3, 9.6 and 20.5 days), the inner 
planet being the lightest planet detected to date around a Sun-like star. Again, 
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Fig. 13. Left: Phase-folded orbits for the three planets around HD 69830. Right: Radial- 
velocity curve as a function of time. Panels a and b show the cumulative signal of the 
three planets. Panel c shows the orbit of the outer planet when the observations are 
binned to one data point per observing run. Figures from ILovU et all l|2006l) . 

one should note the very small semiamplitudes of the three reflex motions (2.0, 
2.5, and 4.6 ms^^) and the low rms of the residuals around the best-fit model 
(0.85 ms-i). 

4.2.2 Nascent statistics and emerging trends 

Occurrence rate of low-mass planets The recent discoveries of planetary 
systems with Neptune-mass planets or super-Earths confirm that low-mass planets 
are common in short-period orbits around solar- type stars. A first estimate based 
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on the HARPS survey suggests that 30 ± 10% of the G-K d warfs harbor a planet 



with TUp sini < 30 Mg, and P < 50 days (jMavor et al.ll2008f) . 



Multiple planet systems with low-mass planets The fraction of multiple 
planet systems with at least one low-mass planet is very high (~80%). One can 
argue that this high value stems from an observational bias, the detection of some 
low-mass planets resulting from a special interest for the host star and its first- 
found planet (s). However, this is less and less true, and many of the recent de- 
tections are the result of systematic surveys, not of follow-up observations. The 
fraction of multiple planet systems with at least one low-mass planet is thus prob- 
ably intrinsically high. 



Mass distribution As shown in Fig. [HI the low-mass planets discovered recently 
seem to build up a new population which distinguishes itself from the population 
of giant planets. That is, the mass distribution looks bimodal and the decrease 
previously observed below 1 Mjup may not be due entire ly to se l ection effects. 
Remarkably, the synthesis population models of iMordasini et al. ( 2007 ) predict 
a planetary initial mass function which possess a local minimum near ^40 M0 
(0.13 Mjup) and a local maximum for Neptune-mass planets (0.05 Mjup). Present 
data on low-mass planets are thus consistent with theoretical expectations and 
future discoveries will put to the test the relative overabundance of Neptune-mass 
planets predicted by theory. 



Planet-metallicity correlation The correlation observed between the occur- 
rence rate of giant planets and stellar metalhcity seems less pronounced for low- 
mass planets ( Udry et al. 20061 : Bonfils et al. 2007 ). If confirmed, this trend could 
indicate that most of the short-period Neptune-mass planets are subcritical plan- 
etary cores rather than the remnants of evaporated giant planets. Future re- 
sults from Doppler surveys and from transit observations (Sect. [5]) should provide 
important constraints to characterize the nature and the origin of short-period 
Neptune-mass planets. 



Stellar mass dependence M dwarfs seem to harbor more low- m ass than giant 
planets in short-period orbits (jEndl et al. l l2008l : lBonfils et al.ll2007^ ■ implying that 
the ratio of Jupiter- to Neptune-mass planets may decrease with decreasing stellar 
mass. The sample of low-mass planets is still too small to draw r obust conclusions , 
but aga in this trend is consistent with theoretical predictions teenz et al. 1 120081: 
Ida fc L in 2005; Laughlin et al. 2004). The planet population orbiting M dwarfs is 
thus probably quite different from the planet population orbiting solar-type stars. 



4.3 Exploring the limitations on Doppler measurements with HARPS 

In the context of the GTO program, significant efforts have been put into charac- 
terizing the limitations on HARPS measurements with the aim of identifying the 
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main obstacles to the detection of Earth twins. The two factors that mainly limit 
the instrumental precision of HARPS are the stability in the illumination of the 
spectrograph and the wavelength calibration. To minimize variations in the illu- 
mination pattern, a specific guiding algorithm has been implemented on the ESO 
3.6-m telescope. This algorithm controls the position of the photocenter at the 
0.05 arcsec level, keeping the guiding noise below 30 cms~^. As to the wavelength 
calibrat ion, it has been imp roved recently and the global uncertainty is now of 20 



(jLovis fc Pepell200<f l. 



Stellar noise is another major limitation on very high precision Doppler mea- 
surements. As shown in Fig. [71 solar-like oscillations are clearly detected with 
HARPS. The strategy adopted to minimize this oscillation noise consists in setting 
the exposure times to 10-15 minutes to average out the signal over 2-3 oscillation 
periods. Although time consuming, this approach allows to keep the residual noise 
below ~50 cms~^ for late-G and K dwarfs. On intermediate and long time scales, 
Doppler measurements are affected by stellar granulation and stellar activity, re- 
spectively. The impact of stellar granulation has not been quantified precisely yet. 
As to stellar activity, present results indicate that well-s elected quiet G-K dwarfs 
exhibit a jitter level below 1 ms~^ ( Pepe fc Lovid 2008f ). Moreover, binning the 



observations over time scales comparable to the rotational periods of the stars per- 
mits to significantly average out stellar noise. Proceeding in this way, a precision 
of 30-40 cms~^ has been obtained on a few test cases (e.g. HD 69830). Further 
work will be needed to quantify the ultimate limitation imposed by stellar activ- 
ity, but a careful selection of chromospherically inactive targets probably allows 
to achieve a precision better than 50 cms^"'^ on many stars. 

The short-term (1 night) precision of HARPS has been characterized by aster- 
oseismology observations of the star a Cen B carried out during the commissioning 
phase. Subtracting photon noise (17 cms^^) and stellar noise (44 cms^^) from 
the global rms of the measurements (51 cms^^) leaves 20 cms~^ for all other er- 
ror sources (guiding errors, instrumental errors, influence of the atmosphere, etc.). 
On short time scales, the precision is thus mainly limited by photon noise (tar- 
gets from planet search programs are fainter than aCenB) and by the intrinsic 
stability of the stars. 

The long-term precision of HARPS is more difficult to characterize. The his- 
togram of the observed radial-velocity dispersions for the stars from t he high- 



preci sion GTO program (G and K dwarfs) presents a mode at 1.4 m s ( M ayor fc Udrv 



20081) . Yet, part of the observed scatter is likely caused by undetected multiple 



planet systems with low-amplitude Doppler signals and by stellar activity. The 
long-term performances of HARPS arc thus probably better illustrated by consid- 
ering the residuals around the best orbital solutions published so far. For several 
stars (e.g. HD 69830, HD 40307), these residuals are below 1 ms~^, which indicates 
that HARPS reaches a sub-ms~^ precision on the long term. If reproducible on a 
sufficiently large sample of stars, the precision of 35 cms"^ obtained for HD 69830 
would show that we are not far from detecting Earth-mass planets in the habitable 
zones of Sun-like stars. 
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4.4 Perspectives and future instruments 

The exciting results obtained with HARPS have motivated new studies to push 
down the limits of Doppler spectroscopy to reach the 1 cms~^ precision level. 
From the instrumental perspective, the experience gained with HARPS indicates 
that reaching this precision level sho uld be possible, though several issues will 



have to be solved ( Pepe fc Lovid 2008 ). From the scientific perspective, the main 
limitation on ultra high Doppler precision will probably be stellar noise. But again 
the results obtained with HARPS are encouraging, and applying an adequate 
observing strategy on carefully selected targets offers good prospects of averaging 
out stellar noise down to 10-20 cms~^ on short and intermediate timescales for 
some dozens of stars. Future ultra-stable high-resolution spectrographs are thus 
under stu dy. The major Euro pean project is CODEX for the Extremely Large 



Telescope ( Pasquini et al.llioOSh . CODEX aims at reaching a precision of 1 cms 



over at least 10 years to directly measure the expansion of the Universe. For planet 
searches, CODEX would provide an unprecedented facility for the radial-velocity 
follow-up of Earth- twin candidates detected through photometric transits (see next 
section). As an intermediate step between HARPS and CODEX, E SO foresee s 



the realization of ESPRESSO for the Very Large Telescope ( D'Odoricd I2007D . 



ESPRESSO would be a second-generation HARPS-type instrument designed to 
achieve a long-term precision of 10 cms~^. These characteristics would allow to 
carry out a systematic search for Earth-mass planets. 



5 Transiting planets 

In the course of their orbit, some extrasolar planets intersect the line of sight 
between the star and the observer, momentarily occulting part of the starlight. 
These regular passages in front of the stellar disk are called transits, while pas- 
sages behind the stellar disk are called occultationsH. For a distant observer, 
the probability to see a transiting geometry is low for long-period planets such 
as Jupiter or the Earth (0.1% and 0.5%, respectively), but is higher for close-in 
planets such as 51Pegb (10%). The prospects of detecting extrasolar planets by 
looking for transits have thus considerably improved with the discovery of hot 
Jupiters. To date, 52 extrasolar planets are known to transit their parent stars. 
These planets are very precious because they provide us with unique information 
about their physical properties and atmospheres. 



5.1 Observations of planetary transits and occultations 
5.1.1 Photometric transits 

Two examples of transit light curves are shown in Fig. 1141 Besides their repetition 
time (which corresponds to the planet's orbital period), such light curves provide 



^Transits are sometimes called primary eclipses, while occultations are also referred to as 
secondary eclipses or antitransits. 



Detection and Characterization of Extrasolar Planets 



31 




^ J -0.04 -0-02 0.02 0-04 

-0.04 -0.02 0.03 0.04 

Phase Phase 



Fig. 14. Phase-folded transit light curves for HD 189733 b obtained from the ground 
(lejf) and from space (right). Features A and B observed in the residuals of the right- 
hand light curve are attributed to the occultation of starspots by the planet. Figures 
from iBouchv et all l|2005ch and lPont et"aD (|2007l ). 



four observables: the transit duration, the transit depth, the ingress/egress dura- 
tion, and the central curvature. These two last parameters correspond to relatively 
subtle features in the curves and can be measured only on high signal-to-noise data. 
By supplementing the transit depth, duration, and the orbital period with some 
external knowledge of the stellar mass and radius, we can estimate the inclina- 
tion of the planetary orbit and the planet's radius (e.g. Seager &: Mallen-OrnelasI 
20031) . However, transit photometry provides no information about the planet's 



mass. Transit photometry and Doppler spectroscopy are thus highly complemen- 
tary in their outcomes, and combined observations allow to determine both the 
true mass and the radius of a planet. 

Transiting planets can be discovered and characterized through two different 
paths. In the first case, the planet detection is made through Doppler measure- 
ments and follow-up photometric observations reveal the transit. The first extra- 
solar planet found to trans i t its host star, HP 209458 b. was disco vered in this way 
(|Charbonneau et alllioool : iHenrv et al]l200d iMazeh et aDl2000l ). In the second 
case, photometric observations detect a transit and follow-up Doppler measure- 
ments confirm the planetary nature of the transiting object. Spectroscopic con- 
firmation is essential in this case bec ause transit photometry suffers from a high 
rate of false alarms (e.g. lBrownll2003[ ). Photometric detection is the most efficient 
channel to discover transiting planets, but the bright and nearby stars of Doppler 
surveys are better suited for complementary studies (Sect. ISTO]) . 

The passage of Jupiter in front of the Sun (or of a Neptune-like planet in 
front of an M5 dwarf) produces a dip of 1% in the light curve, while a tran- 
sit of the Earth produces a dip of 0.01%. Photometric transits from short- 
period Jupiter-sized planets can be detected from the ground, but transits from 
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Earth-sized planets are detectable only from space (see e.g. Brown 20081 for 
further details). Numerous photometric surveys searching for planetary tran- 
sits are presently active and the rate of discovery has increased rapidly over the 
past two years. Amo ng grou nd-based surveys, H AT (iBakos et al.ll2007t). OGLE 



(|Udalski et a l.l'2002al). TrES (lAlonso et al.ll2004l) . WASP (|Cameron et al.l[2007h . 
and XO (jMcCu lloug h et al.l 120061 ) have successfully detected transiting planets. 
These surveys typically monitor tens of thousands of stars with a photometric 
precision better than 10 mmag. CoRoT - the first space-borne mission largely 
dedicated to the discovery of pla netary transits - h as also found several planets 
since its launch in December 2006 ( Barge et al. 2008f) . During its 2.5-year mission, 
CoRoT will monitor ^^60,000 stars for 150 continuous days, permitting the detec- 
tion of super-E arths (>2 Rq) and larger planets with periods below '^50 days (e.g. 
Moutoul[2006h . 



An example of photometric follow-up: GJ 436 The low-mass planets dis- 
covered since 2004 have prompted intensified follow-up photometry to check for 
possible transits. The first - and currently only - Neptune-sized planet found 
to transit its parent star is GJ 436 b, who se transits were de tected at a non- 
professional observatory in the Swiss Alps ( Gillon et al. l2007bl ). This discovery 



opened a new window for planetary studies a nd triggered immedia te follow-up 
observations with the Spitzer Space Telescope (jGillon et al ][2007a ; Deming et al 



[20071 : iDemorv et al.H2od7i see also Sect. 15X31) . The mass and radius of G J 436b 
suggest an internal structure similar to that of Neptune, with a rocky core sur- 
rounded by a thick water layer and a thin gaseous envelope (Fig. [16]). Nonetheless, 
without additional observationa l constraints, the in ternal composition of GJ 436 b 
cannot be determined in detail ( Adams et al. 20081 see also Sect. I5.2.2p . 



An example of Doppler follow-up: the OGLE candidates In 2002-2003, 
OGLE was t he first photometric survey to deliver a series of planet candidates 
( Udalski et a l. 2002a, c,b, 2003), and Doppler follow-up programs were immedi- 
ately initiated to identify genuine planets. The task was difficult, however, be- 
cause most of the OGLE targets are too faint {V G [14;18]) to be observed with 
standard planet-search spectrographs. The Geneva team resorted to using the 
FLAMES/UVES multi-fiber spectrograph on the Very Larg e Telescope (VLT) to 



follow u p spectroscopically many of the OGLE candidates ([Bouchv et al. I l2005bl : 



IPont et a l. 2005). This effort has been very successful, leading to the discov- 
ery of 3 of the 5 transiting planets fo und among the 137 first OGLE candidates 
([Bouchv et al ] \2004 [Font et al.[[2004[ ). The follow -up campaign has also yielded 
empirical occurrence rates for the main types of false alarms: transits by small 
stars (38%), multiple star systems (24%), false transit detections (14%), and graz- 
ing eclipses (7%). Flanets represent only ^--^7% of the total number of candidates. 

The follow-up of the OGLE candidates has demonstrated the necessity and 
the efficiency of Doppler observations in identifying genuine planets among the 
numerous false alarms. Yet, characterizing the transits exhibited by the faintest 
OGLE targets is close to being impractical. The OGLE survey is thus near the 
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Fig. 15. Phase-folded radial-velocity curve for HD 147506 showing the Rossiter- 
McLaughlin effect (inset). Figure from iLoeillet et alJ l|2008l ). 

upper limiting magnitude for spectroscopic follow-up with existing ground-based 
facilities. 

5.1.2 Spectroscopic transits: the Rossiter-McLaughlin effect 

A planetary transit produces not only a photometric signal but also a spectro- 
scopic signal called the Rossiter-McLaughlin effect. This anomaly in Doppler shift 
measurements traces the passage of the planet in front of the rotating stellar disk 
(Fig. fTS)) . That is, when the planet occults part of the approaching (blueshifted) 
stellar hemisphere, the integrated starlight appears slightly redshifted, an effect 
interpreted as an anomalously high radial velocity. The reverse effect is observed 
when the planet occults part of the receding stellar hemisphere. For extraso- 
lar planets, the Rossiter-McLa ughlin effect was first observed and modeled for 
HD 209458 ('Oueloz et al. Since then it has been detected for 11 transiting 

systems (see .Win n 2QD8 for a review). Two important parameters can be obtained 
by modeling the Rossiter-McLaughlin effect: the direction of planetary revolution 
(prograde or retrograde orbit), and the angle A between the sky projections of the 
orbital and the stellar rotation axes. This angle gives us some information about 
the degree of alignment of the planet's orbital angular momentum vector with 
the stellar spin axis (spin-orbit alignment), which provides in turn an important 
constraint for the formation models of hot Jupiters (Sect. 

5.1.3 Precise monitoring of transits and occultations 

Close-in planets with large planet-to-star radius ratios and transiting nearby bright 
stars permit a variety of follow-up studies which provide us with deeper insight 
into the physical properties of these planets. The two most favorable systems 
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for detailed studie s are HD 209458 b dHenrv et al.ll2000HCharbonneau et"alll2000l: 
Mazeh et al.ll2000l ) and HD 189733 b (|Bouchv et al.ll2005cl ). 

Transit light curves of extreme accuracy - as obtained with the Hubble Space 
Telescope, e.g. right panel of Fig. [14] - reduce the ambiguity between the es- 
timated parameters (stellar radius, planetary radius, orbital inclination, stellar 



Brown et al. 2001 



limb darkening) and yield very precise planetary radii (e.g 
[Pont et al.l l2007h . High-quality transit time series allow additionally to search 
for evidence of planetary satellite s, Saturn-type rings, or additional p lanets (e.g. 



Sartoretti fc Schneideil[l99 9: Bar nes k Fortnevll2004t lAgol et al.ir2005l ) 



Follow-up observations using the HubbleEl and SpitzeiQ Space Telescopes have 
provided the first data on the atmospheres of extrasolar planets. At visible and 
ultraviolet wavelengths, the light emitted by the planet is negligible and the tech- 
nique of transmission spectroscopy has been employed to study the chemical coin - 
position and the physical properties of planetary atmospheres (e.g. BrownllioOll ). 
Transit spectroscopy uses the fact that the major sources of opacity in a planetary 
atmosphere depend on wavelength. Due to this dependence, the planet's effective 
radius - and thus the transit depth, which is the ratio of the planet-to-star surface 
area - vary with wavelength. Searching for relative changes in transit depth as 
a function of wavelength is thus a means to probe the absorption properties of a 
planet's atmosphere (e.g. Charbonneau et all 20021 : Vidal-Madiar et al. 20031 ). 

At infrared wavelengths, the planet-to-star contrast ratio improves consider- 
ably, permitting a direct detection of the photons emitted from the planet. Since 
hot Jupiters are believed to rotate synchronously with their orbital period (i.e. 
they are tidally locked), they present permanent day and night sides. The depth 
of the secondary eclipse (occultatio n) allows to characterize the day-side flux (ther - 
mal emission) from the planet fe.g. lDeming et aL l l2005HCharbonneau et al.1l2005l) . 
while observations at various orbita l phases allow to characterize the planet's lon - 



20061 : iKnutson et all l2007t ) 



gitudinal temperature profile (e.g. [Harrington et al 
In addition, the timing and duration of the secondary ecl i pse place useful con - 
straints on the orbital eccentricity (e.g. iDeming et al.l 120071 : iDemorv et al.ll2007l ). 
Up t o now, the thermal emis sion from six extrasolar planets has been detected 
(e.g. ICharbonneau et al.ll2008l . and references therein). 



5.2 Noteworthy results obtained from transiting planets 
5.2.1 A new class of planets: the very hot Jupiters 

Three of the five planets di scovered among the first 137 OGLE can didates have 
orbital periods of 1.5 days ( Konacki et al. 2003t Bouchv et al ] l2004h . which origi- 
nally placed them in an empty region of the period distribution, below the pile-up 
of hot Jupiters. These new planets were accordingly termed very hot Jupiters. 
Since the sensitivity of Doppler spectroscopy increases with decreasing period, the 



^http: / /bubble. nasa.gov 
^http:/ /www. spitzer.caltech.edu 
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absence of very hot Jupiters in the discoveries of radial-velocity surveys looked 
surprising. I G audi et al. showed that transit surveys are actually ~6 times 

more sensitive to planets with periods of 1 day than of 3 days, a fact that could 
explain the observations. Since then, Doppler surveys have found several planets 
with periods of 2-3 days, while transit surveys have detected many regular hot 
Jupiters. These discoveries have bridged the gap between the two populations and 
confirm that there is no significant incompatibility between the results of Doppler 
and transit surveys. 



5.2.2 Mass-radius relationships and bulk compositions 

The measurement of a planet's mass and radius give us the mean density, from 
which we can infer some information about the planet's bulk composition (Fig.ll6p. 
Most of the known transiting planets have densities below 1.5 gcm^"^, so they 
must be composed primarily of hydrogen and helium, like Jupiter and Saturn. 
Historically, this has been an important confirmation of their planetary status. 

Some transiting giant planets have high densities which cannot be explained 
by a pure H/He composition. The most extreme case is HD 149026b (indicated 
on Fig. [T6|l . whose small radius requires t he presence of ^ 70 M^ of heavy ele- 
ments (2/3 of the total mass) in its interior ( Sato et al. 2005f) . Giant planets with 
small radii are accounted for due to variations in the amount of heavy elements 
sequestered in a central core and/q r distributed throughout the gaseous envelope 
(|Guillotll20dilBurrows et al.ll2007a[ ). Since the presence of a solid core is a natural 
characteristic of formation by core accretion, giant planets with small radii appear 
to support the core accretion model. Nonetheless, if they are common, planets 
with 2/3 of their mass of heavy elem ents are a challenge to any formation model, 
inclu ding core accretion (jSato et al.i i2005t llkoma et all l2006t iBroeg fc Wuchterll 
2003). 

Contrary to HD 149026 b, some giant planets such as HD 209458 b (also in- 
dicated on Fig. I16[) have unexpectedly large radii and low densities. Accord- 
ing to standard models of irradiated planets, these hot Jupiters are too big f or 
their mass and age, a characteristic which is still unexplained (see lGuillotll2008l or 
Burrows et al. 2007al for further details). 

Low-mass planets are expected to display a wide variety of compositions and 
internal structures . According to theoretical predictions fe.g. iFortnev et al.ll2007 : 
Seager et al.ll2007h . Neptune-mass planets can be composed primarily of refrac- 
tory materials (rock, iron, planetary ice sEI, but also of H/He, or of a mixture of 
these possibilities. As to super-Earths, they are presumably composed of refrac- 
tory species. It should be noted that there is some ambiguity and degeneracy in 
inferring bulk compositi ons from mean dens ities, especially for planets with a mass 
between 5 and 20 M^ ( Adams et al. 2008^. For example, mass and radius mea- 



surements alone are insufficient to distinguish between an ocean planet (a planet 



^Planetary ices comprise ammonia, methane, and water (in any phase, not necessarily solid). 
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Fig. 16. Mass-radius diagram for planets. Transiting hot Jupiters (diamonds), planets 
from the Solar System (dots) , and GJ 436 b are indicated. Figure from iGillon et al.l 
l|2007bl '). adapted from the models of iFortnev et al.1 |20o3). 



with >25% of water by mass) and a rocky planet with a massive H/He gaseous 
envelope. 



5.2.3 Spin-orbit alignment 

Rossiter-McLaughlin observations obtained so far show that in all except one case 
transiting p lanets are on prograde orbits aligned with the stellar rotation axis (see 
Winnll2008l for a review) . These results suggest that the migration of hot Jupiters 
generally preserves spin-orbit alignment, an observation that disfavors scattering 
processes or Kozai migration as the main formation channel for hot Jupiters. 
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5.2.4 Planetary atmospheres 



Hot Jupiters are tidally locked and receive >10,000 times more radiation from 
their stars than Jupiter does from the Sun. Due to this high and asymmetric irra- 
diation, their atmospheric properties are expected to differ significantly from those 
of Jupiter and Saturn. Atmosphere models of hot Jupiters predict transmission 
spectra characterized by strong absorption lines due to the sodium and potassium 
resonance doublets at visible wavelengths, and by stro ng molecular bands of water, 
carbon monoxide, and methane in the infrared fe.g. iBrownl 120011 ). A significant 
source of uncertainty in atmosphere models is the possible presence of hazes or 
clouds, which can block the stellar flux at any height in the atmosphere. Another 
key question regarding hot Jupiters is whether atmospheric circulation efficiently 
redistributes the absorbed stellar radiation from the day side to the night side. 

Follow-up studies on HD 209458 b have revealed the presence of so dium, hydro- 

gen, water, and more tentatively oxygen and carbon in its atmosphere (ICharbonneau et al 



2002Hyidal-Madiar et al.ll2003ll2004lBarmanl2007HBurrows et al.l2007bHKnutson et all 

Hydrogen absorption extends beyond the planet's Roche lob^, indicating 
that HD 209458 b has an extended and evaporating atmosphere (see the contribu- 
tion by Ehrenreich) . The infrared broadband emission spectrum of HD 209458 b 
differs significantly from the predictions of standard atmosphere mo dels, which was 



interpreted as evidence for a n atmospheric temperature inversion (jKnutson et al 



2008at iBurrows et"al]l2007bh . The day-night contrast of HD 209458 b h as not been 



measured directly, but atmo spheric circulation seems relatively efhcient (jCowan et al 



2007t iKnutson et al.l[2008bh . 



tection of water in its atmosohere (.Tinetti et al.l 


2007: Ehrenreich et al. 


2007 


Beaulieu et al. 2008; Swain et all 20081 


Charbonneau et al.l 2008t Barman 


2008 


Fortnev & MarlevI 2007; Grillmair et al. 


2007). The good agreement obtained re- 



cently between the broad band emi s sion s pectrum measured by Charbonne au et al 



( 2008f l and the models of BarmanI (2008) provides convincing evidence that water 
absorption does shape the infrared spectru m of HD 18 9733 b. Methane and sodium 
have also been detected in its atmosphere ( Swain et al . 2008; Redfield et al. 200§)- 
The emergent spectrum of HD 189733 b is well matched by standard models whic h 
do not include an atmospheric temperature inversion (^ Charbonneau et al.ll200"8h . 
HD 189733 b is only slightly hotter on the day side than on the night side, indi- 
cating that the energy fro m the irradiated hemisphere is efficiently redistributed 
throughout its atmosphere (*Knutson e t al. 2007 . 2008a ). 

Detailed follow-up studies on HD 209458 b and HD 189733 b have thus revealed 
both similarities and fundamental differences in their atmospheric properties. To- 
gether with the data obtained on a few other transiting systems (see Burrows et al.|[2008 
for a review), the above results support the idea that the atmospheres of hot 
Jupiters diverge into two groups: those with and those without a temperature 



^The Roche lobe outlines the volume surrounding an object within which material is gravita- 
tionally bound to it. 
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inversion (jFortnev et al.ll2008l : iBurrows et al. I l2008h . 



5.3 Perspectives and future instruments 

In 2009, CoRoT will be joined by the Keplei0 satellite, which wiU monitor 100,000 
main-se quence stars for 4 continuous years to detect the transits of Earth- like plan- 
ets (e.g. Basri et al. 2005[) . Besides the discovery of hundreds of low-mass planets. 



CoRoT and Kepler will provide the first measure of the occurrence of Earth- and 
Neptune-sized planets at short and intermediate periods. The frequency of Earth- 
sized planets orbiting in the habitable zones of solar-type stars is a key parameter 
for the design of the next generation of project missions, which will aim at char- 
acterizing the atmospheres of Earth-like planets. 

Space-based transit searches do not eliminate the need for spectroscopic follow- 
up, but the higher quality of the photometry reduces the rate of false alarms. 
The radial-velocity follow-up of the CoRoT planet candidates is conducted with 
CORALIE, SOPHIE, HARPS, the Coude spectrograph on the Alfred- Jensch tele- 
scope, and FLAMES/UVES for the faintest targets. The spectroscopic follow-up 
of the Kepler planet candidates will be much more challenging and will require a 
larger investment of telescope time. To perform the Doppler follow-up of interest- 
ing planet candidates identified by Kepler, the Harvard University and the Geneva 
Observatory have joined in a collaboration to build and operate the HARPS-N 
spectrog raplllll- This improved HARPS-type instrument will be installed on the 
William-Herschel Telescope at La Palma Observatory in 2009. In synergy with 
Kepler, HARPS-N is a promising step towards the detection of rocky planets in 
the habitable zones of solar-type stars. 

Regarding the precise monitoring of transits and occultations, Spitzer will run 
out of cryogen in spring 2009, which will eliminate its ability to observe at wave- 
lengths longer than 4.5 /im. An extended "warm" Spitzer mission has been pro- 
posed, in which the spacecraft would continue to operate (at the same sensitivity) 
in the two shortest wavelength channels. Unlike Spitzer, Hubble should be reborn 
with the servicing mission scheduled for early 2009. Beyond Hubble and Spitzer, 
we will have to wait the launch of the James Webb Space Telescoptfll in the mid- 
dle of the next decade to continue and extend the study of extrasolar planetary 
atmospheres. 

6 Conclusion and perspectives 

For centuries, our knowledge of planetary systems has been based on a single ex- 
ample, the Solar System. During the last thirteen years, Doppler spectroscopy has 
allowed the detection of ~300 extrasolar planets, opening a new window on plan- 
etary sciences. These discoveries have considerably broadened our appreciation 



"http:/ /kepler. nasa.gov 

-"^http: / /obswww.unige.ch/Instrumcnts/harps_north 
^http:/ /www.j wst.nasa.gov 
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of the diversity of possible planetary systems and have revolutionized our ideas 
about planet formation. 

By the standards of the Solar System, many extrasolar planets exhibit singu- 
lar properties, including very short orbital periods, high eccentricities, and large 
masses. Stars of different masses seem to harbor different planet populations and 
multiple planet systems appear to be the rule. For solar-type stars, stellar metal- 
licity is a strong predictor of the presence of giant planets in the inner regions of 
the system. Migration via tidal interaction with the gaseous protoplanetary disk 
appears to play a central role in the formation and early evolution of these giant 
planets. Some extrasolar planets were foimd in environments relatively hostile to 
their presence (pulsars, spectroscopic binaries), indicating that planet formation is 
a robust process. The recent discoveries of Neptune-mass planets and super-Earths 
at close (Dopplcr spectroscopy) and intermediate (microlcnsing) separations sup- 
port the theoretical prediction that low-mass planets are more numerous than their 
giant cousins. 

Transiting planets are only a small subset of all extrasolar planetary systems, 
but they offer an unique opportunity to deepen our understanding of planet struc- 
ture and composition. With the exception of one hot Neptune, transiting planets 
are gas giants with a wide range of masses and radii. Follow-up observations using 
the Hubble and Spitzer Space Telescopes provided the first glimpse into the at- 
mospheres of some hot Jupiters, and the first detections of photons emitted from 
extrasolar planets. Planetary emission spectra reveal fundamental differences in 
atmospheric properties, leading to the suggestion that the irradiated atmospheres 
of hot Jupiters bifurcate into two groups. 

The above observations provide an ever-growing database with which planet 
formation theories can be put to the test. By combining observed properties 
of protoplanetary disks, models of planet formation and evolution, and observa- 
tional selection effects, theorists now generate in a self-consistent way synthetic 
planet populations which can be compared with the actual detections. Recent 
models based on the core accretion paradigm successfully reproduce key observa- 
tional features such as the mass distribution and the planet-metallicity correlation. 
Although much remains to be understood about planet formation, these new the- 
oretical achievements reinforce the idea that most of the known planets formed 
through core accretion. 

To progress in our understanding of planet formation and to appreciate the 
significance of our Earth, we must detect a wide range of planetary systems, in- 
cluding planets and host stars both similar and dissimilar to our own. This can be 
achieved through two different paths: by further developing existing techniques to 
increase their capabilities, and by developing new techniques to explore different 
regions of the parameter space. The characterization of transiting planets illus- 
trates the benefit of combining the results obtained with two different observational 
techniques. 

As we have seen, Doppler spectroscopy has reached neither its instrumental 
nor its astrophysical limits. The latter will probably be dominant, but applying 
an adequate strategy on carefully selected targets offers good prospects of averag- 



40 



Title : will be set by the publisher 



ing out stellar noise down to 10-20 cms ^ on short and intermediate time scales 
for some dozens of stars. The next generation of HARPS-type instruments aim- 
ing at a Doppler precision of 10-20 cms'^ (HARPS-N, the ESPRESSO project) 
should therefore be able to detect Earth-mass planets in short-period orbits, and 
to characterize some of these planets up to periods of ~1 year. Ultra-stable spec- 
trographs aiming at a Doppler precision of 1 cms~^ should be technically feasible. 
If stellar noise permits, such instruments would allow the detection and the precise 
characterization of Earth twins. 

While Doppler spectroscopy is intrinsically sensitive to short-period planets, 
astrometry and direct imaging are particularly suited to detect planets in wide 
orbits. Up to now, these two techniques have very marginally contributed to plan- 
etary studies, but this is going to change thanks to the upcoming deployment 
of new facilities. In 2009, the PRIMAI^'^I instrument for the VLT interferome- 
ter should become fully operational, allowing the astrometric de tection of planets 



down to the mass of Neptune in the separation range 1-5 AU (jLaunhardt et al 



2008h . From space, the GAl43 mission (scheduled for launch in 2011) will detect 



thousands of extrasolar planets, revolutionizin g the database o f extrasolar plan- 
ets with known orbital elements and masses ( Casertano et al.l [2008). Although 



astrometric observations are sufficient to determine three-dimensional orbits and 
true planetary masses, additional Doppler measurements will be of great help in 
constraining the large number of free parameters, especially for multiple planet 
systems. Regarding direct imaging, a few "extreme" adaptive optics sy stems are 



under construction, including SPHERE for the VLT (jBeuzit et al.ll2008l ) and GPI 



for the Gemini Observatory ( Macintosh et al. 2008[ ). Starting in ~2011, these 



instruments will detect, characterize the atmospheres, and determine the radii 
of wide (^--^5-100 AU) and massive (>l-5 Mj^p) planets of various ages ('--^lO-lOOO 
Myr) . All together, these upcoming facilities will probe the intermediate and outer 
regions of protoplanetary disks, providing new information about the formation, 
evolution, and migration of giant planets. 

Doppler spectroscopy and microlensing have opened the way to the detection 
of low-mass planets, and mature detection techniques now aim at identifying an 
Earth-mass planet located in the habitable zone of its host star. Microlensing 
surveys can yield the frequency of such plane ts, but they do not a llow a precise 
characterization of their planet candidates (e.g. Beaulieu et aLllioOsh . Space-based 



transit photometry is presently the most sensitive method for planetary detection, 
and the pending Kepler mission offers good prospects of finding the first Earth- 
sized planet potentially habitable. Yet, Doppler spectroscopy is a serious competi- 
tor, and it may win the race. In any case, both techniques will play a key role in 
detecting and characterizing Earth-mass planets in the near future. Due to their 
ubiquity and low mass, M dwarfs are ideal targets for the detection of terrestrial 
planets. In addition to optical radial velocities, infrared Doppler spectroscopy rep- 
resents an attractive new technique to search for and characterize planets around 



^■^http: / /obswww.unige.ch/Instruments/PRIMA 

^^http:/ /gaia.esa.int/science-e/www/area/index.cfm?fareaid=26 
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late-type (>M4) M dwarfs (e.g. iMartm et a"l]|2005l or the SPIRou projec10). In 



the more distant future, space-borne astrometry could contribute to the detection 
of low- mass planet s as well. The SIM Litcflfl project would serve this purpose (e.g. 



Tanner et al. but its future is highly uncertain. 



If it is approved, the warm Spizter mission will allow the community to retain 
the ability to characterize the atmospheres of some transiting planets for a few 
more years. In the middle of the next decade, JWST will take over and will allow 
to study the atmospheres of Earth-sized planets located in the habitable zones of 
main-sequence stars. On a longer time scale, infrared space interferometers could 
be built to directly image potentially habitable Earth-like planets and to search 
for life signatures in their atmospheres. 

In summary, Doppler spectroscopy will continue to serve as a front-line de- 
tection technique for the next years. It will additionally provide an increasingly 
important support to alternative planet search techniques, in particular transit 
photometry and astrometry. The golden age of radial velocities is thus far from 
its end and Doppler spectroscopy promises additional exciting results to come. 
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